Chaotic Zeeman effect: a fractional diffusion-like approch

It is shown that the chaotic Zeeman effect of a quantum system can be formally viewed as a result of fractional calculus. The fractional calculation brings into the equations the angle \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta $$\end{document}θ formed between the internal and the external magnetic field applied to the atom. The further the fractional coefficient \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}α is from the ordinary case corresponding to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =1$$\end{document}α=1, the more important the chaotic effect is. The case corresponding to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =1$$\end{document}α=1 does not depend on the angle \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta $$\end{document}θ, obtaining the nonchaotic situation known in the literature. Non-Gaussian distributions correspond to non-stationary variables. Considering a Lorenzian type distribution, we can make a connection between the fractional formalism and random matrix theory. The connection validates the link between fractional calculus and chaos, and at the same time due to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta $$\end{document}θ angle, it gives the phenomenon a physical interpretation.

Fractional equations model anomalous diffusion in many systems.In the paper 1 , it is shown that the probability density function of tracer particles' radial displacements is strongly non-Gaussian showing algebraic decaying tails, and with the help of the fractional diffusion model it can be reproduced the shape and space-time scaling of the probability.The paper 2 studied front dynamics in reaction-diffusion systems where anomalous diffusion is due to asymmetric Levy flights.The approach consists of replacing the Laplacian diffusion operator by a fractional diffusion operator.Motivated by work on contaminant spreading in geological formations, in the paper 3 , the authors propose and investigate a fractional advection-diffusion equation describing the biased spreading packet.In the solution of the variable-order fractional diffusion equation, the authors of the paper 4 identified a new advection term that causes ultraslow spatial aggregation of subdiffusive particles due to dominance over the standard advection and diffusion terms in the long-time limit.In 5 , the authors derived a fractional Fokker-Planck equation for subdiffusion in a general space-and time-dependent force field from power law waiting time continuous time random walks biased by Boltzmann weights.The experiment does not involve quantum mechanics, it involves liquid crystal dynamics, and thus the related but still unsolved problem of whether it is indeed a quantum problem.
The Zeeman effect is the consequence of placing an atom in a uniform external magnetic field, and results in the splitting of a spectral line into several components.The nature of the Zeeman splitting depends critically on the strength of the external magnetic field compared to the internal one.When the external magnetic field B is small compared to the internal one B int , the fine structure dominates and the problem is treated perturbatively.In the situation B ≫ B int the Zeeman effect is dominant and the fine structure becomes perturbation.Also called the normal Zeeman effect, this phenomenon was first explained by Lorentz using Bohr's atomic model.However, when the forces due to the external magnetic field become comparable to the Coulombian binding forces acting on the electron, the known theories break down and the energy levels become irregular, as we represented in Fig. 1.This transformation of the atomic spectrum from regular to irregular coincides with the situation in which the corresponding classical system transits from the regular to the chaotic regime.Next, we title this image of well-stirred spaghetti, chaotic 6 .
In the treatment of the non-chaotic Zeeman effect, we use the fact that due to the strength of the external magnetic field, the internal magnetic field is parallel to it, and the resulting theory does not depend on the angle θ that these two magnetic fields make between them.Inspired by the fact that the sequence of nearest level spacings has similarities with the diffusion process of a particle 7 , in this work, we propose an approach based on fractional calculus.Any deviation of the fractional coefficient α from α = 1 , makes the theory depend on the angle θ .In this paper, we show that this angle is responsible for the chaotic behavior of the process.We also note that random matrix theory predicts that the probability distribution of the distance between energy levels provides an excellent description of the energy levels in atoms excited in strong magnetic fields.We show that this is not accidental, by simply identifying the mean level density function with the Lorenzian profile we obtain an oscillating behavior similar to that given by the fractional calculus.
The Riemann-Liouville fractional integral is the inverse operator of the fractional derivative in the Caputo sense.Chaos is understood as an extreme sensitivity of the solutions of a system of differential equations to the initial conditions.Because Caputo derivatives share this sensitivity, they are a useful tool for extending regular behavior to irregular behavior.It is also known that chaos cannot appear in continuous systems with a number of equations less than three, except for fractional-order systems 8 .The Caputo derivative is also successfully used in optimal control problems to remove chaotic behaviors 9 .Last but not least, it should be mentioned that the Caputo derivative has the property of preserving the history of interactions in the dynamic behaviors it simulates, a characteristic property of singular kernel derivatives 10 .
The paper is structured as follows: in Section "Theoretical investigation" we conduct the theoretical investigation of the Zeeman effect by changing the definition of mechanical work from the Riemann integral to the Riemann-Liouville fractional integral.This allows us to link the energy levels of the atoms in the magnetic field to the θ angle.In Section "Results analysis", we analyze the new behavior of the Zeeman effect and make the connection between fractional calculus and random matrix theory.We conclude the paper with the observation that the theory presented in this article can be successfully extended to any type of chaos, regardless of whether it is classical or quantum.

Theoretical investigation
Definition 1 11 Magnetic moment of the atom in a uniform magnetic field µ l is a consequence of the fact that the electron performs an orbital movement around the nucleus, and is defined as where L is the orbital angular momentum of the electron, e is the electron charge and m is the mass of the electron.
Definition 2 11 When an atom is placed in a magnetic field B , the torque acting on it is defined It means that the magnetic field works on the dipole to rotate it.Definition 3 11 The mechanical work done by the field on the system will be stored in terms of potential energy, and it must correspond to the rotation of the dipole from π/2 to θ Definition 4 We make the following ansatz where the Riemann-Liouville integral operator I α is defined as 12 Proposition 1 13 We have the following result where the two-parameter Mittag-Leffler function is given by the following expression 14 (1)  3) and using Eq. ( 2), we obtain Because L z = m l , we get the desired result.
Remark 1 When α = 1 , F(θ, α = 1) = 1 , regardless of θ .This means that in the ordinary case, the potential energy does not depend on the angle between L and B.
Proof We start from the fact that we get and as a consequence In Fig. 2, we represented the function F(θ, α) , for different values of α and θ .It can be seen that if α = 1 , the function is not dependent on θ and is equal to 1. Definition 5 11 In the absence of the magnetic field, the jump of an electron between two energy levels E 1 and E 2 is accompanied by a photon with frequency where E 2 is the excited energy.Definition 6 If the atom is placed in a magnetic field then the total energy E t of the atom is Proposition 3 For an atom in a magnetic field, the transition of an electron is accompanied by a photon of frequency with m l = m 2 l − m 1 l .
Proof According to Eq. ( 4), the energies corresponding to the two energy levels involved are B F(θ, α), This means that in the situation where m l = 0 , chaos has no contribution.Remark 3 When α = 1 , we have F(θ, α) = 1 and U(θ) = µ B m l B , a result that allows us to obtain the ordinary situation from the literature 11

Results analysis
We are interested in analyzing the behavior of Eq. ( 4), more precisely the dependence of U α (θ) according to the magnetic field B and the function F(θ, α) , knowing that m l = 1 and that µ B = 5.788 × 10 −5 eV T −1 .
The angle θ will oscillate around 0, and for a value description we assume that we have the following vari- ation of θ with the time θ(t) = sin 5t .In Fig. 3 we represented U α (θ)/µ B , and as we can see, by how much α departs from the value 1, the oscillation depending on θ becomes stronger.So the value of α is closely related to the deviation of the energy spectrum from the ordinary behavior.The farther α has a value from 1, the more prominent the chaotic behavior.
Fixing the value of α = 0.8 , In Fig. 4 we analyze the behavior of U α (θ)/µ B depending on the magnetic field B and the angle θ considering the same time dependence as in Fig. 2.This tells us that for a given value of B, the function can have any of the values given by t ∈ [0, 1.5] s.In this way, a spectrum is obtained that shows irregular patterns that resemble well stirred spaghetti 6 .
Random matrix theory provides an excelent description of the spacing of the measured energy levels in excited nuclei.Although there are no random elements in atomic phenomena, the statistics of the energy-level spacing agreed remarkably well with the prediction of the random matrix theory.The unfolded enegy levels are defined 7 where ρ(β, E ′ ) is the mean level density function, and β stands for some parameters defining the functional form of ρ.
We notice that Eq. ( 4) comes from an integral over time, and Eq. ( 6) is described by an integral over energy.If there is a connection between the two functions, this is given by the Fourier transform.
Non-Gaussian distributed data arise when the mean statistics vary with time, which means that the variable could be non-stationary: time-varying (its trend changes with time), or cross-correlated (it changes depending B F(θ, α). on its previous value or the values of other variables).Since the fractional calculation takes into account the history of the process, we think of using a non-Gaussian distribution.We choose a Lorentzian distribution for the unfolded energy levels, i.e., and then its Fourier transform has the expression and the real part of this function when Ŵ → 0 is a sinusoidal.This behavior is similar to that represented in Fig. 3, for a given value of α .The fact that a link can be created between the two formalisms once again justifies the usefulness of fractional calculus in describing chaotic behaviors.

Conclusions and outlook
Inspired by the description of diffusion, we used fractional calculus to describe the chaotic Zeeman effect.Fractional calculus has recently been used to model various phenomena.In the paper 15 , the authors studyed the existence of numerical solution and stability of a chemostat model under fractal-fractional order derivative.The analysis proved to be an excellent tool for discussing the fractal characteristics of porous media.The manuscript 16 focused on the modeling and numerical solution of the dynamical model of typhoid fever.The modeling was based on the Atangana-Baleanu operator with the Mittag-Leffler function in Caputo sense.In the article 17 , the authors studyed the Drinfeld-Sokolov-Wilson equation considered in fractal-fractional sense with exponential  www.nature.com/scientificreports/decay and Mittag-Leffler type kernels.The dynamics of COVID-19 is investigated in 18 , where the authors used a fractional order SEIR model.Compared to random matrix theory, in which we do not have random elements in atomic physics, in our theory the chaos depends on the fluctuation of the θ angle between the intrinsic magnetic moment and the mag- netic field applied to the atom.Figure 3 shows us that when α = 1 , case corresponding to the non-chaotic Zee- man effect, the energy levels do not depend on the θ angle.Eq. ( 5) represents the chaotic unfolding of the energy spectrum.As can be seen in Fig. 3, the fluctuation amplitudes of the energy levels are modulated by the fractional value of the coefficient α .In Fig. 4, you can see the irregular spectrum reminiscent of well stirred spaghetti.
Since the fractional calculation introduces the history of the process into the dynamics, we choose to use a non-Gaussian distribution for the data.The Fourier transform of the Lorentz distribution has a sinusoidal shape, which makes the connection with the oscillating behavior described in Fig. 3.The link between the two formalisms justifies the effectiveness of fractional calculus in describing chaotic processes.At the same time, due to the physical interpretation of the angle θ , the formalism described by us gets a physical interpretation.
We expect the oscillation frequency of L around B to be high, or in other words θ to oscillate rapidly around the value θ = 1 , value corresponding to the non-chaotic energy spectrum.This property is known as rigidity of the energy spectrum.This effect is similar to the vibrations obtained in high voltage wires.Analog 7 , we can consider the energy spectrum as a discrete signal and the sequence of energy levels as a time series.In fact, we proceeded similarly in Fig. 3, where we generated a time series when we considered θ(t) = sin 5t .In a time series, antipersistence means that an increasing or decreasing trend in the past makes the opposite trend in the future, and quantum systems are very antipersistent (characterized by 1/f noise).The rigidity of the spectrum as well as its rapid oscillation around the non-chaotic value obtained with our theory, are consistent with the antipersistent cataloging from 7 .
Turbulence is easily associated with fractals, and fractals represent the geometry of chaos.Probably the best known example of turbulence is the vortex on Jupiter.Looking closer and closer, we see first larger vortices, then smaller ones, keeping the structure of the whole on a reduced scale.Geometric objects with this structure are called fractals and have a different geometry than ordinary geometry.For example, a sphere viewed more and more closely becomes more and more flat.Fractals are different, when we enlarge them we see more and more structure.
An idea to study the Zeeman effect from this point of view is to modify Eq. ( 3), replacing the Riemann-Liouville integral operator with an fractal or fractal-fractional integral operator.In the literature there are several operators of this type 19,20 .Maybe this type of operators bring out a new dynamic of the electronic spectrum highlighting aspects still undiscovered related to the Zeeman effect.
In fact, the analysis in this paper, although applied to quantum chaos, should work successfully in any kind of chaos, regardless of whether it is classical or quantum.

Figure 1 .Proposition 2
Figure 1.Chaotic distribution of energy levels: the spectra show irregular patterns that resemble well-stirred spaghetti.